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\^ , We report the discovery of a planetary system from observation of the high-magnification microlensing event 

■ OGLE-2012-BLG-0026. The lensing light curve exhibits a complex central perturbation with multiple features. 

■ We find that the perturbation was produced by two planets located near the Einstein ring of the planet host star. 
We identify 4 possible solutions resulting from the well-known close/wide degeneracy. By measuring both the 
lens parallax and the Einstein radius, we estimate the physical parameters of the planetary system. According 
to the best-fit model, the two planet masses are 0.11 Mj and 0.68 Mj and they are orbiting a G-type main 
sequence star with a mass ^ 0.82 Mq. The projected separations of the individual planets are beyond the snow 
line in all four solutions, being ^ 3.8 AU and 4.6 AU in the best-fit solution. The deprojected separations 
are both individually larger and possibly reversed in order This is the second multi-planet system with both 

k> , planets beyond the snow line discovered by microlensing. This is the only such a system (other than the Solar 

. System) with measured planet masses without sin / degeneracy. The planetary system is located at a distance 

■ 4.1 kpc from the Earth toward the Galactic center It is very likely that extra light from stars other than the 
lensed star comes from the lens itself. If this is correct, it will be possible to obtain detailed information about 
the planet-host star from follow-up observation. 
Subject headings: planetary systems - gravitational lensing: micro 
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ABSTRACT 



1. INTRODUCTION 

According to the standard theory of planet formation, rocky 
planets form in the inner part of the protoplanetary disk of 
a star, within the snow line, where the temperature is high 
enough to prevent th e condensation of water and other sub- 
stances onto grains (Ravmond et al. 2007). This results in co- 
agulation of purely rocky grains and later in the formation of 
rocky planets. On the other hand, formation of giant planets 
starts at some distance from the host star, beyond the snow 
line, where the low temperature of the protoplanetary disk al- 
lows for condensation of water ice, enabling more rapid ac- 



cumulation of solid material into a large planetary core. The 
core accretes surrounding gas and later evolves into a gas gi- 
ant planet (Pollack et al. 1996). Due to the difference in the 
conditions of formation, rocky terrestrial and gas giant plan- 
ets fo rm at different regions of planetary systems (llda & LinI 
12004)). 

Planet formation theories were initially based entirely on 
the Solar System. Planet search surveys using the radial ve- 
locity and transit methods have discovered numerous planets, 
including many gas giants. Most of the giant planets have 
short-period orbits and are thus much closer to their parent 
stars than the giant planets in the Solar System. These dis- 
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coveries have inspired new theoretical work, such as vari- 
ous migration theories to resolve the contradiction between 
the present-day and assumed-birth locations of giant planets. 
However, it is essential to discover more planets located be- 
yond the snow line to bridge the gap between the Solar System 
and exo-planet systems. 

Microlensing is an important tool that can provide observa- 
tional evidence to check planet formation theories. The plan- 
etary signal in a microlensing event is a short-term perturba- 
tion to the smooth standard light curve of the lensing event 
induced by the planet's host star (Mao & Paczyriski 1991; 
iGould & Loeblll992D . The perturbation is caused by the ap- 
proach of a lensed star (source) close to the planet-induced 
caustics. The caustic is an optical term representing the en- 
velope of concentrated light that forms when light ray is re- 
fracted by a curved surface, e.g. a curved region of bright 
light when light shines on a drinking glass filled with wa- 
ter or a rippling pattern at the bottom of a swimming pool. 
For gravitational bending of light, caustics represents the po- 
sitions on the source plane at which the lensing magnifica- 
tion of a point source becomes infinity due to the singularity 
in the lens-mapping equation. When a caustic is formed by 
a planetary system, it forms multiple sets of closed curves 
each of which is composed of concave curves that meet at 
cusps. Among them, one set is always located close to the 
planet's host star (central caustic) and other sets are located 
away from the star (planetary caustic). Due to these caustic 
locations, a perturbation induced by a central caustic (cen- 
tral perturbation) occurs near the peak of a high-magnification 
event, while a perturbation induced by a planetary caustic oc- 
curs on the wings of a lensing light curve. The caustic size, 
and thus the planet detection sensitivity, is maximized when 
the planet is located in the region near the Einstein ring of the 
host star, which is often referred to as the lensing zone. For 
a typical Galactic lensing event, the physical size of the Ein- 
stein radius is re ^ 3.Q{M /Mq)^I'^ AU, where M is the mass 
of the lens. This is sim ilar to the snow Une of 2.1(M /Mp,y 
AU, w here 0.7 < p < 2 dlda & Linl2005l: [Kennedy & KenvonI 
12008'). Therefore, microlensing is an efficient method for de- 
tecting cool planets located at or beyond the snow line. 

Microlensing is also sensitive to multiple -planet systems. 
Microlensing detection of multiple planets is possible because 
all planets in the lensing zone affect the magnificat ion pattern 
around the primary lens ( Griest & Safiza deh'199S). and thus 
the signatures of the multiple planets can be de tected in th e 
light curve of a high magnification event (IGaudi et al.lll999^ ■ 
Indeed, the first two-planet microlensing syste r n was discov- 
ered through this channel (iGaudi et al.l 120081 : iBennett et alJ 
l2OT0h . 

In this paper, we report the second discovery of a multiple- 
planetary system composed of two planets located beyond the 
snow line discovered by using the microlensing method. 

2. OBSERVATION 

OGLE-2012-BLG-0026 was discovered by the Optical 
Gravitational Lensing Experiment Early Warning System 
(OGLE EWS; Udalski 2003) on 2012-Feb-13 toward the 
Galactic bulge at (a,6)mm = (17''34'"18^.70,-27°08'33".9), 
i.e., il,b) = (0.195°, 3.066°), based on /-band observations 
with the 1.3m Warsaw telescope at Las Companas Obser- 
vatory in Chile. At HJD' (= HJD-2450000) = 5989.3 (Mar 
2), the Microlensing Follow-Up Network (/iFUN) issued an 
alert, saying that the event would reach high magnification 
(A > 100) 2.5 days hence, and urging immediate observations. 



Given this extremely early date (3.5 months before the target 
reaches opposition), and hence very short observing window 
from any individual site, this quick alert was crucial to mo- 
bilizing near-24-hour coverage from multiple sights. /^FUN 
covered the peak with observations from 6 telescopes: the 0.4 
m at Auckland Observatory, 0.4 m at Farm Cove Observatory 
(FCO), 0.4 m Possum Observatory, 0.4 m Turitea Observa- 
tory in New Zealand, 0.4 m Klein Karoo Observatory (KKO) 
in South Africa, and 0.3 m Perth Exoplanet Survey Telescope 
(PEST) in Australia. The Auckland data were taken in R band, 
and the remainder were unfiltered. In addition, ptFUN ob- 
tained supplementary V /I/H color observations using the 1 .3 
m SMARTS telescope at Cerro Tololo Inter- American Obser- 
vatory (CTIO). Follow-up observation was stopped after the 
perturbation, but normal monitoring of the event was contin- 
ued by the OGLE survey until the event returned to its base- 
line brightness, and beyond. 

Photometric reductions of data were carried out using pho- 
tometry codes developed by the individual groups. The OGLE 
data were reduced by a photometry pipeline that applies the 
image subtraction method based on the Difference Image 
Analysis technique ( Alard & LuDton 1998) and further devel- 
oped by Wozniak (2000) and Udalski (2003). Initial reduction 
of the /iFUN data wer e processed using a pi peline based on 
the DoPHOT software (ISchechter et al.lll993l). The data were 
reprocessed by using the pySIS package (lAlbrow et al.ll2009l) 
to refine the photometry. 

Figure m shows the fight curve of OGLE-2012-BLG-0026. 
The magnification of the lensed source star flux reached 
Amax ^ 129 at the peak. The perturbation is localized near 
the peak region of the light curve. It exhibits a complex struc- 
ture with multiple features. At a glance, the perturbation is 
composed of two parts. One part, at HJD' ~ 5990.2 has a 
positive deviation with respect to the unperturbed single-lens 
light curve. The other part, centered at HJD' ~ 5991.2 has 
a negative deviation. Such a complex perturbation pattern is 
unusual for central perturbations. 

3. INTERPRETATION OF LIGHT CURVE 

To interpret the observed light curve, we initially conducted 
modeling based on a standard two-mass lens model. Descrip- 
tion of a binary-lens light curve requires 3 single lensing pa- 
rameters plus another 3 parameters related to the binary na- 
ture of the lens. The single lensing parameters include the 
moment of the closest lens-source approach, to, the projected 
lens-source separation at that moment, mq, and the Einstein 
time scale, fg, of the event. The Einstein time scale is defined 
as the time required for the source to cross the angular Ein- 
stein radius of the lens, 9e- The binary lensing parameters 
are the projected binary separation in units of the Einstein ra- 
dius, 5, the mass ratio between the lens components, q, and 
the angle between the source trajectory and the binary axis, 
a. Since central perturbations usually involve caustic cross- 
ings or approaches, an additional parameter of the source ra- 
dius in units of 6e, (normalized source radius), is needed to 
account for the attenuation of lensing magnifications by the 
finite size of the source. A central perturbation in the light 
curve of a high magnification event can be produced either 
by a planetary companion positioned in the lensing zone or 
by a binary companion located away from the Einstein radius 
dHan & Gaudil2008l;lHanl2009b . From a thorough search con- 
sidering both possible cases of central perturbations, however, 
we could not find any solution that explains the observed per- 
turbation. 
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Fig. 1 . — Light curve of OGLE-2012-BLG-0026. The lower panel shows the whole view of the light curve and the upper panel shows an enlarged view of the 
perturbation. Data points taken from different observatories are marked by different colors along with the passbands of observation. The passband notation "N" 
represents that no filter is used. The dotted and solid curves represent the solutions from single and triple lensing modeling, respectively. The two residuals in the 
lower panel are presented to show the difference between the fits of the models with and without the parallax effect. 



We theref ore conside r the possibility of an additional lens 
companion. lHan| ( l2005h pointed out that many cases of cen- 
tral perturbations induced by multiple planets can be approx- 
imated by the superposition of the single-planet perturba- 
tions in which the individual planet-primary pairs act as in- 
dependent binary lens systems. Based on this "binary su- 
perposition" approximation, we then conduct another binary 
lens modeling of the two light curves each of which is pre- 
pared by including only a single perturbed region. The first 
light curve includes the negative perturbation region during 
5990.5 < HJD' < 5992.0 and the second includes the positive 
perturbation region during 5989.0 < HJD' < 5990.5. From 
this, we find that the two partial light curves are well fitted by 
binary lensing light curves with two different planetary com- 
panions. 

Based on the parameters of the individual planetary com- 
panions, we then conduct triple-lens modeling, which requires 
to include three additional parameters related to the second 
companion. These are the normalized projected separation 
from the primary, s^, the mass ratio between the second com- 
panion and the primary, q^, and the position angle of the sec- 
ond planet with respect to the line connecting the primary and 



the first planet, ■\\}. We denote the parameters related to the first 
planet as s\ and q\. We note that the first and second planets 
are designated as the ones producing the negative and posi- 
tive perturbations, respectively. With the initial values of the 
triple-lens parameters obtained from the binary superposition 
approximation, we locate the minimum using the Markov 
Chain Monte Carlo (MCMC) method. 

We find that consideration of the parallax effect is needed 
to precisely describe the light curve. The parallax effect is 
caused by the change of the observer's position due to the 
orbital motion of the Earth around the Sun during the event 
(Gould 1992). For OGLE-2012-BLG-0026, the event time 
scale is of order 90 days, which comprises a significant frac- 
tion of the Earth's orbital period, and thus parallax can affect 
the lensing light curve. We therefore introduce two additional 
parallax parameters of tte.w and tte.b, which represent the two 
components of the lens parallax vector tte projected on the 
sky in the north and east equatorial coordinates, respectively. 
The direction of the parallax vector corresponds to the lens- 
source relative motion in the frame of the Earth. The size of 
the parallax vector is defined as the ratio of Earth's orbit to 
the physical Einstein radius projected on the observer plane. 
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Fig. 2. — Four possible geometries of the lens systems. Each of the left 
panels shows the relative positions of the two planets with respect to the host 
star. Coordinates are centered at the position of the host star and all lengths 
are normahzed by the Einstein radius. The circle with a radius 1.0 represents 
the Einstein ring. Also marked is the trajectory of the source star (a curve with 
an arrow). Right panels show an enlarged view of the central region around 
the host star. The cuspy feature represents the caustics. The two small empty 
dots represent the source star at the moments of the major perturbations. The 
size of the dots is scaled to represent the source size. 

i.e. tte = A\]/[0e/(DI^-D^^)], where Dl and Ds are the dis- 
tances to the lens and source star, respectively. We find that 
the parallax effect improves the fit by Ax^ ^ 240. We note 
that "terrestrial" parallax effect due to n on-cospatial observa- 
tories on the Earth (Gou ld et alj 12009b is not detected since 
the separation between the 2 observatories covering the caus- 
tic crossings (Auckland and PEST) is not big enough. 

For high-magnification events, planets described by {s,q) 
and {s~\q) induce nea rly identical caustic s and thus nearly 
identical light curves jGriest & Safizadehl |1998). This is 
known as the close/wide degeneracy in lensing light curves. 
We, therefore, investigate other solutions with initial values 
of the projected separations based on the first solution. From 
this search, we find 4 sets of solutions with {si < 1, ^2 < 1), 
(si > 1, S2 < 1), (si < 1, 52 > 1), and (^i > 1, 52 > !)■ We 
denote the individual models as "A", "B", "C", and "D". 

If the source were exactly on the ecliptic, it would exactly 
obey the "ecli ptic degeneracy", wh ich takes (uo,a,TTE.±) — > 
-{uQ,a,'!TE,±) (ISkowron et al.ll20lll) . Here tte.x is the compo- 
nent of tte perpendicular to the ecliptic at the projected posi- 
tion of the Sun at fo, which in the present case is almost ex- 
actly north. OGLE-2012-BLG-0026 lies < 4° [(A,/3)y2ooo = 
(264. 28°, -3. 83°)] from the ecliptic. Hence it is not surpris- 
ing that we find almost perfect ecliptic degeneracy for all four 
solutions, with Ax^ < 1- We do not show these redundant 
solutions. 

In principle, one can detect lens orbital motion, and even 
if not, tangential or bital motion can be degenerate with tte,± 
dBatista et al.ll201 fl) . which would degrade the parallax mea- 
surement. We fit for all possible orientations of co-planar cir- 
cular orbital motion (hence two additional parameters to de- 



fine the orbit plane). However, we find, first that there is no 
improvement in the fit (Ax^ ^ 1)> and second that there is no 
degeneracy with tte x. We therefore do not further consider 
orbital motion. 

In Table |2l we present the four solutions. In Figure |2l we 
also present the geometries of the lens system corresponding 
to the individual solutions. We find that both planets respon- 
sible for the perturbation are located near the Einstein ring 
of the host star. The negative perturbation was produced by 
the passage of the source trajectory on the backside of the 
central caustic induced by a planet located almost on the Ein- 
stein ring of the host star, while the positive perturbation was 
produced when the source passed the sharp front tip of the 
arrowhead-shaped central caustic induced by another planet 
located at a position slightly away from the Einstein ring. The 
position angle of the second planet with respect to the first 
planet is ip ^ 220°. We note that although the close/wide de- 
generacy causes ambiguity in the star-planet separation, the 
mass ratios of the two degenerate solutions are similar to each 
other We find that the mass ratios are ^1 ^ L3 x lO"'* and 
7.8 X lO""*. 
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4. PHYSICAL PARAMETERS 



The combined measurements of the angular Einstein ra- 
dius and the lens parallax enable an estimate of the lens mass 
and distance. With these values, the mass and distance are 
determined by M = 6e/{k,tte) and the De = AU/(7rE6'E + 7rs), 
where n = 4G/( c^AU) and tt s = AU/Ds is the parallax of 
the source star (Gould 2000). By modeling the Galactic 
bulge as a bar, the distance to the source star is estimated by 
Ds = DGc(cos^+sin^cos 4)/ sinc^), where D pc = 8.16 pc is the 
Galactocentric distance ( Nataf et al.ll2012l) . (. = 0.195° is the 
Galactic longitude of the source star, = 40° is the bar orien- 
tation angle with respect to the line of sight. This results in 
Ds = 8.13 kpc, corresponding to tts = 0.123 mas. 

The Einstein radius is given by 9e = where 6** is the 

angular radius of the source, which we evaluate from its de- 
reddened color (y-/)o and brightness /q u sing an instrumental 
color-magnitude diagram dYoo et al.ll2004)) . We first measure 
the offset A[(y-/),/] = (-0.07,2.79) mag of the source from 
the clump, whose de-reddened position is known indepen- 
dently, r(V-/),/ 1n. = (1.06, 14.45) mag (iBensbv et al.ll2011t 
iNataf et al.l2012l) . This yields the de-reddened color and mag- 
nitude of the source star [(y-/),/]o,s = (0.99 17.23) mag 
Thus t he source is a bulge subgiant. Note that iBensbv et al.l 
(120131) estimate a very similar color (V -/)o = 1 .00 mag based 
on a high-resolution VLT spectrum taken at high magnifica- 
tion. Figure |3] shows the locations of the source star and the 
centroid of the giant clump in the color-magnitude diagram 
that is constructed based on the OGLE data. Once the de- 
reddened V-/ color of the source star is measured, it is trans- 
la ted into Y —K co lor by using the V -I versus Y -K relations 
of Bes sell & Bretd (ll988) and then the angular source radius 
is estimated by using the rel ation between the Y - K color and 
the angular radius given bv lKervella eTaTI (120041) . The mea- 
sured angular source radius is 0^, = 1.55 ±0.13 /xas. Then, the 
Einstein radius is estimated as Qe = 0.91 ±0.09 mas. 

In Table [U we list the physical parameters of the planetary 
system determined based on model D, which provides the best 
fit. Other models are similar. According to the estimated mass 
of ^ 0.82 Mq, the host star of the planets is a G-type main 
sequence. The individual planets have masses of ~ 0.11 Mj 
and 0.68 Mj. Their projected separations are ~ 3.8 AU and 
4.6 AU, respectively, and thus both planets are located well 



HAN ET AL. 



5 



14 



18 



20 



~~| — I — I — I — I — I — I — I — !~ 

giant clump 
source 

blend 



1.5 



2.5 



3.5 



Fig. 3. — Locations of the source star and blend in the color-magnitude 
diagram. Also marked is the position of the giant clump centroid that is used 
to calibrate the color and magnitude of the lensed star. Black dots are stars in 
the field where the source star is located. 

beyond the snow line 1.8 - 2.4 AU of the host star. The plan- 
etary system is located at a distance ^4.1 kpc from the Earth 
toward the Galactic bulge direction. 

To be noted is that the blend is likely to be the lens it- 
self. Blended light refers to the extra light coming from stars 
other than the lensed star The source of blended light can 
be unresolved stars located close to the lensed star or the 
lens itself. The fraction of blended light among total ob- 
served light is measured from modeling of the light curve. 
Considering the distance to the lens ^4.1 kpc, one can as- 
sume that the lens and the source experience roughly the same 
amount of reddening and extinction. Under this approxima- 
tion, the estimated de-reddened color and brightness of the 
blend determined with the reference of the giant clump are 
(y-/,/)o.b = (0.93, 16.62) mag. Then, under the assumption 
that the blend is the lens combined with the distance modulus 
^ 13.05 mag to the lens, the estimated absolute magnitude 
of the blend is M/b ^3.6 mag. Then, the color and the ab- 
solute magnitude of the blend roughly agree with the color 



(V — I)o ^ 0.9 mag and the absolute magnitude Mj ^ 4.2 mag 
of a G5V star, that corresponds to the lens, suggesting that 
blended light comes from the host star Moreover, from dif- 
ference imaging the magnified and unmagnified images, we 
find that the blend and source have identical positions within 
measurement error 50 mas), further tightening the identi- 
fication of the blend with the lens. 

This identification opens many possibilities for future 
follow-up observations. AO imaging could confirm and re- 
fine the host star flux measurement by isolating it from fainter 
ambient stars. At / = 19 mag, the host is bright enough to 
obtain a spectroscopic metallicity measurement, which would 
be the first for a microlensing planet host. Finally, late-time 
astrometry could improve the mass and distance determina- 
tions by refining the 9e measurement from the magnitude of 
the proper-motion vector ( Alcock et al. 2001 ) and refining the 
tte measurement from its direction (iGhosh et al.ll2004l) . 

TABLE 1 
Physical Parameters 



parameter 


quantity 


mass of the host star 


0.82 ± 0.13 M0 


mass of the first planet 


0.11 ±0.02Mj 


mass of the second planet 


0.68±0.10Mj 


projected separation to the first planet 


3.82 ± 0.30 AU 


projected separation to the second planet 


4.63 ± 0.37 AU 


distance to the planetary system 


4.08 ± 0.30 kpc 
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MICROLENSING MULTIPLE-PLANET SYSTEM 



TABLE 2 
Lensing Parameters 



parajnetcrs 


model A 


model D 


model C 


model U 


X^/dof 


2689.7/2667 


2681.0/2667 


2680.8/2667 


2674.4/2667 


(0 (HJD') 


5991.52 ±0.01 


5991.52 ±0.01 


5991.52 ±0.01 


5991.52 ±0.01 


Mo 


-0.0094 ± 0.0001 


-0.0097 ± 0.0001 


-0.0094 ± 0.0001 


0.0092 ± 0.0001 


tE (days) 


92.45 ± 1.08 


89.21 ±0.17 


92.85 ± 1.11 


93.92 ± 0.58 


■«l 


0.957 ± 0.001 


1.034 ± 0.001 


0.957 ± 0.001 


1.034 ± 0.001 


qi (10-^) 


1.32 ±0.02 


1.37 ±0.01 


1.30 ±0.02 


1.30 ±0.01 


■«2 


0.812 ±0.004 


0.819 ±0.004 


1.257 ±0.006 


1.254 ±0.006 


?2 (10-^) 


7.97 ± 0.26 


7.84 ± 0.25 


8.04 ± 0.25 


7.84 ± 0.21 


a 


1.283 ±0.001 


1.284 ±0.001 


1.284 ±0.001 


4.999 ± 0.001 


v> 


2.389 ± 0.002 


2.387 ± 0.001 


2.391 ± 0.001 


3.891 ± 0.001 


(10-3) 


1.73 ±0.02 


1.81 ±0.01 


1.72 ±0.02 


1.72 ±0.01 




-0.029 ± 0.054 


-0.097 ± 0.012 


0.001 ± 0.028 


-0.072 ±0.051 




0.123 ±0.009 


0.137 ±0.003 


0.123 ± 0.005 


0.114 ±0.004 



Note. — HID'=HJD-2450000. 



